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I . INTRODUCTION 


The  use  of  fire  as  a weapon  by  man  precedes  recorded  history. 
Anthropologists  theorize  that  early  man  used  burning  brush  to  scare 
wild  animals  from  his  habitat.  More  recently,  the  ancient  Romans  and 
Greeks  used  firebrands  and  boiling  oil  or  water  to  repel  invaders. 

During  the  early  twentieth  century,  many  incendiary  and  pyrotechnic 
munitions  were  developed.  Three  decades  ago,  the  world  was  introduced 
to  the  atomic  bomb  and  its  sun-like  fireball.  Given  the  resourceful- 
ness of  man  and  the  knowledge  of  modem  science,  the  "ray  gun"  may  be- 
come reality. 

All  of  these  weapons  have  the  same  basic  purpose:  the  incapacitation 
and/or  death  of  the  enemy  by  pair,  and  burning.  In  the  past,  weapons 
were  fabricated  out  of  locally  available  materials  that  were  obtained 
at  minimal  cost.  Because  the  monetary  investment  was  so  low,  a trial- 
and-error  approach  was  used  in  developing  such  new  weapons  as  were  re- 
quired. During  the  past  fifty  years,  however,  weapons  and  weapons 
systems  have  become  both  omplex  and  expensive.  Because  of  this,  man 
has  turned  to  science  both  to  develop  the  weapons  and  to  guarantee 
their  effectiveness  when  developed.  The  most  expeditious  method  of 
determining  a weapon's  effectiveness  is  to  test  a prototype  model  against 
the  intended  target.  However,  the  intended  target  in  these  cases  is 
man  and  the  indiscriminate  burning  of  humans  is  proscribed  by  most 
cultures.  There  exist  several  alternatives:  testing  on  human  volunteers 
under  carefully  controlled  conditions,  testing  on  human  surrogates 
(animals),  and  mathematical  simulation.  The  first  alternative  is  used 
to  some  extent,  however,  the  temperature  range  over  which  testing  is 
done  is  limited  by  the  reluctance  of  the  volunteers  to  undergo  permanent 
bodily  damage.  The  use  of  animals,  a general  practice,  suffers  from 
several  disadvantages:  (1)  the  physiology  of  the  animal  is  not  identical 
to  that  of  the  human,  (2)  animals  are  expensive  to  buy,  feed,  and  house, 
and  (3)  the  administrative  procedures  required  to  gain  approval  for 
animal  experiments  are  formidable.  The  last  alternative,  mathematical 
simulation,  has  been  used  to  a limited  extent  during  the  past  two  dec- 
ades, The  equations  developed  were  strictly  empirical  relationships 
that  attempted  to  predict  subcutaneous  skin  temperatures  on  the  basis 
of  incident  irradiance.  Because  of  the  complexity  of  the  subject,  many 
simplifying  assumptions  were  made.  The  model  described  in  this  report 
views  the  subject  as  a simple  heat  transfer  problem  amenable  to  analysis 
by  theoretical  means. 
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II.  MODELING  TECHNIQUES 


In  order  to  model  any  system  successfully,  it  is  necessary  to 
establish  a relationship  between  the  response  of  the  system  and  one 
or  more  parameters  that  can  be  determined  either  empirically  or  theoret- 
ically. Thus,  the  question  to  be  answered  for  this  model  was  "Is  there 
a relationship  between  skin  temperatures  and  skin  bums?".  This  question 
was  answered  by  Hardy1  who  made  extensive  experimental  observations  of 
bum  thresholds  as  a function  of  skin  temperature.  In  order  to  calculate 
skin  temperatures,  however,  a knowledge  of  the  thermophysical  properties 
of  the  skin  and  underlying  tissue  is  needed.  NASA2,  as  part  of  the 
space  program,  has  compiled  extensive  data  on  human  responses  to  the 
aerospace  environment.  The  needed  thermophysi cal  data  was  found  in 
these  volumes. 

Figure  1 presents  a cross  sectional  view  cf  a portion  of  the 
volume  of  flesh  described  by  this  model.  The  skin  varies  from  0.5mm 
on  the  eyelids  to  5.0mm  on  the  back.  Since  the  tissue  underlying  the 
skin  serves  as  either  a heat  source  or  a heat  sink,  a volume  2.5cm  in 
depth  has  been  modeled.  The  epidermis,  which  is  approximately  500um 
thick,  is  avascular  and  does  not  significantly  affect  the  occurrence  of 
first  and  second  degree  bums.  Therefore,  the  temperature  of  the  outer- 
most layer  of  the  dermis  at  a depth  of  SOOnm  was  selected  as  the  basis 
for  the  bum  criteria. 

The  bum  criteria,  taken  from  Hardy,  are  shown  in  Figure  2.  The 
area  to  the  upper  right  of  the  solid  line  denotes  irreversible  tissue 
damage, i.e.,  a second  degree  bum;  the  area  between  the  dashed  and 
solid  lines  denotes  transient  erythema,  i.e.,  a first  degree  bum. 
Histologically,  a first  degree  bum  is  characterized  by  immediate  pain, 
which  continues  after  exposure  and  by  ensuing  redness  (caused  by  dilation 
of  arterioles)  in  the  exposed  area  which  disappears  in  about  24  hours. 

The  first  degree  bum  is  said  to  be  reversible  tissue  injury.  The 
second  degree  bum  is  characterized  by  pain  with  either  no  visible 
effect  or  by  an  immediate  appearance  of  blanching,  loss  of  elasticity, 
swelling  or  blisters.  After  6 to  24  hours,  an  eschar  forms  over  the 
injured  area  and  is  flexible  and  tan  or  brown  if  the  injury  is  light, 
or  thick,  stiff  and  dark  if  the  injury  is  severe.  A second  degree 
bum  does  not  involve  the  full  thickness  of  the  skin  and  the  remaining 
cells  are  able  to  regenerate  a normal  skin.  A third  degree  bum  is 
defined  as  the  destruction  of  the  full  skin  thickness.  All  of  the 
cells  arc  destroyed  and  arc  unable  to  regenerate  a normal  skin. 

A finite  difference  technique  has  been  employed  in  this  model  to 
perform  the  heat,  transfer  calculations.  These  techniques  are  particularly 
suited  to  problems  involving  complex  geometries  and  boundary  conditions, 
non  linearities  such  as  radiation,  and  variable  material  properties. 

The  reader  interested  in  additional  information  on  finite  differencing 
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Figure  1.  Cross  Section  of  Human  Skin 
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is  referred  to  either  Gaskf  or  Meyers4.  The  tissue  to  be  modeled  was 
subdivided  into  the  72  small  volumes  shown  in  Figure  3.  Each  volume  is 
assumed  to  have  uniform  thermophysical  properties  and  temperatures 
throughout.  For  computer  purposes,  a nodal  point,  located  in  the  geo- 
metric center  of  the  volume,  is  assigned  these  properties.  To  effect 
heat  transfer,  mathematical  conductors  are  constructed  to  connect  these 
nodes.  They  are  assigned  the  heat  conduction  properties  of  the  media 
which  they  represent.  The  skin  density  (1.1  g/cm^)  and  heat  capacity 
(0.8  cal/g-°c)  are  taken  from  Reference  3.  The  thermal  conductivity 
data  shown  in  Figure  4 was  taken  from  Weaver5. 

At  this  point  it  is  appropriate  to  discuss  energy  flow  through  tne 
model.  There  are  two  separate  sources  of  thermal  energy  acting  upon 
the  skin:  the  outward  flow  of  heat  from  the  interior  of  the  body  and 
the  inward  flow  of  heat  from  external  sources  incident  upon  the  skin 
surface.  The  first  source  is  fairly  simple  to  model  since  the  rates 
of  heat  loss  for  various  bodily  areas  are  known  from  aerospace  cooling 
studies6.  The  method  used  in  this  case  is  to  input  the  energy  into 
the  innermost  nodes  and  let  the  model  transport  it  to  the  surface  where 
it  is  dissipated  via  convection  and  radiation.  The  second  source  is 
more  complex  since  external  heat  sources  can  be  of  three  types:  radiative, 
convective  and  conductive.  Of  these,  only  radiative  sources  have  been 
considered  to  date.  Subroutines  have  been  written  to  generate  the 
following  pulse  shapes:  square,  delta,  nuclear,  and  laser.  An  additional 
subroutine  allows  the  user  to  input  an  arbitrary  pulse  shape. 

Radiant  sources  have  a spectral  emittance  that  depends  upon  the 
temperature  of  the  source  and  whether  the  source  emits  as  a blackbody 
or  a gray  surface.  This  spectral  distribution  is  important  because 
the  absorptance  of  radiant  thermal  energy  by  human  skin  is  a function 
of  photon  energy.  Subroutines  are  included  that  correct  for  this 
energy  dependence  provided  that:  (a)  a blackbody  temperature  is  known, 

(b)  the  source  spectrum  is  known,  (c)  the  overall  skin  absorptance  is 
known,  or  (d)  in  the  case  of  lasers,  the  emitted  wavelength  is  known. 

If  none  of  the  preceding  are  known,  a worst-case  approach  is  taken 
where  all  energy  is  assumed  to  be  absorbed.  Gray  surface  emission  and 
atmospheric  attenuation  are  not  considered  by  the  model. 


III.  MODEL  VERIFICATION  AND  APPLICATION 


To  date  only  a minimum  of  model  verifications  has  been  performed. 
This  is  not  considered  critical  at  this  time  for  the  following  reasons: 
(1)  the  model  is  still  under  development,  (2)  the  finite-difference 
technique  used  has  been  well  proven,  and  (5)  adjustment  of  constants 
to  correct  for  minor  inaccuracies  is  particularly  easy  in  this  model. 

One  verification  attempt  has  been  made  and  is  shown  in  Figure  5.  The 
peak  calculated  temperature  varied  from  the  peak  experimental  temperature 
by  only  2%.  It  is  felt  that  any  agreement  better  than  this  would  be 
fortuitious. 
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Figure  3.  Nodal  Model  for  Skin  Burn  Calculations. 
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Upon  completion  of  the  model,  additional  verifications  will  be 
performed  using  data  from  carefully  controlled  experiments. 

Future  work  on  the  model  will  include  the  coding  of  conductive 
and  convective  heat  source  routines  to  permit  the  evaluation  of 
pyrotechnic  and  incendiary  weapons  that  bum  by  either  direct  contact 
or  by  the  heating  of  the  ambient  air. 

Because  man  is  the  controlling  factor  in  the  vulnerability 
assessment  of  many  Army  systems,  it  seems  logical  that  the  next  step 
after  this  model  would  be  a vulnerability  code  that  would  predict  the 
degree  of  incapacitation  of  man,  using  data  generated  by  this  model. 
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